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-4— » ' Abstract. We review the theory for photon-photon scattering in vacuum, and some of the proposals 

for its experimental search, including the results of our recent works on the subject. We then describe 
a very simple and sensitive proposal of an experiment and discuss how it can be used at the present 
\ (HERCULES) and the future (ELI) ultrahigh power laser facilities either to find the first evidence 

. ^ . of photon-photon scattering in vacuum, or to significantly improve the current experimental limits. 
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Among other interesting phenomena Oi 0], Quantum Electrodynamics (QED) and 
non-standard theories of the electromagnetic fieldi,ll predict that the vacuum should 
00 \ behave like a kind of virtual electron-positron plasma, thus allowing for photon-photon 

^ ■ (77) scattering. However, the latter prediction has not yet been confirmed by experiment, 

. not even indirectly. In fact, two possible strategies can be used for the quest of 77 
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scattering. On one hand, the cross section for the process will be maximum at a possible 



O ' future photon-photon collider||5|], based on an electron laser producing two beams of 



photons in the MeV range (i.e., having wavelengths in the range of a few fm). A second 
approach will be to perform experiments using ultrahigh power optical lasers, such as 
those that are already in use todaylS fk 

or will be available in the near future[6, |8|], 
^ ■ in such a way that the high density of photons will compensate the smallness of the 

cross section. In this case, the photon energies are well below the electron rest energy, 
and the effect of photon-photon collisions due to the interchange of virtual electron- 
positron pairs can be expressed in terms of the effective Euler-Heisenberg nonlinear 
Lagrangiani^ [I^. This modifies Maxwell's equations for the average values of the 



electromagnetic quantum fields[l 1] and affects the properties of the QED vacuum [|12|] . 

Here, we will briefly review the theory for 77 scattering, and discuss some of the 
proposals for its experimental search. We will then consider the present and future 
ultrahigh power laser facilities, and show that, besides their important technological 
uses, they will provide a unique opport unity fo r the search of 77 scattering, taking into 



account the results of our recent works [|13l.ll4l1. Finally, we will describe a very simple 



and sensitive proposal of an experiment and discuss how it can be used at the present 



(HERCULES) (f] and the future (ELI) [S] ultrahigh power laser facilities either to find 
77 scattering, or to improve significantly the current experimental limits. 



THE QED NONLINEAR WAVE EQUATION 

Even in the absence of matter, photons can interact with each other by interchanging 
virtual electron-positron pairs in a loop. In this sense, the vacuum behaves like some 
kind of virtual electron-positron plasma. We will assume that the photon energy is well 
below the threshold for the production of electron-positron pairs, InigC^ ~ IMeV. This 
means that we will only consider radiation of wavelengths A ^ 2 x 10^ ^^m, which is 
always the case in optical experiments. In this case, the QED effects can be described 
by the Euler-Heisenberg effective Lagrangian density liol fll]. 

^ = ^0 + 5^0 + ^t£y(E ■ B)2, (1) 

being 

^0 = ?fE2-cV) (2) 



2 

the linear Lagrangian density and £o and c the dielectric constant and the speed of light 
in vacuum, respectively. Here, E and B are the expectation values of the electromagnetic 
fields, that can be measured by the usual classical means e.g. using test charges. 

The additional, nonlinear terms, that appear multiplying the parameters and in 
Eq. (HI), are the only two Lorentz-co variant terms that can be formed with the electro- 
magnetic fields at the lowest order above ^q. Therefore, they appear as the first correc- 
tion to the linear evolution in QED in all the alternative theories, such as the Born-Infeld 
theories UM- In QED, ^P""" = J^,^^^ = J^, being 
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This quantity has dimensions of the inverse of an energy density. This means that 
significant changes with respect to linear propagation can be expected when the elec- 
tromagnetic energy density is such that i^p is not very small. While such intensities 
may have an astrophysical or cosmological importance, they are not achievable in the 
laboratory, as we shall discuss in the next section. 

Once the electromagnetic fields are expressed in terms of the four-component gauge 
field = (A*^, A) as B = V A A and E = —cVA^ — the equations of motion are given 
by the Variational Principle: 

0, (4) 
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where F = / ^d^x is the QED effective action. Instead of stud ying th e resulting equa- 
tions for the fields E and B, that can be found in the literature! 111. 1 15|l . it is more con- 



venient to consider the equations for the gauge field components A^. As we discussed 
in Ref. [|T3l.[l4ll. such equations admit solutions in the form of linearly polarized waves. 



e.g. in the x direction, with A" = and A = (A, 0, 0) , provided that: i) the field A does not 
depend on the variable x (a transversality condition) and ii) A{t,y,z) satisfies the single 
equation [fBlfljl: 



= d^A + <9M - d}A + 5eoc2 { (5) 
[[dtAf-3[dyAf-{dAf]d^A+ 
[{dtAY-{dyAY-'i{dAY]d^A- 
[3{d,AY-{dyAf-{d,AY]d^A+ 
4 {d^AdjAd^dtA - d^AdyAd^dyA + dyAdtAdydjA) } , 

where dy = ^„d, = j; and dr = Jj. 

Hereafter, we will restrict our discussion to this case of linearly-polarized solu- 
tion. Note that the plane-wave solutions of the linear Maxwell equations, such as 
^/ cos(k-r— (Ot), where ^/ is a constant, k = (O,^^,,^^) and CO = c|k|, are still solu- 
tions of Eq. (|5]). However, we expect that the non-linear terms proportional to ^, due 
to the QED correction, will spoil the superposition principle. As we shall see, this will 
imply a change in the propagation of crossing waves. On the other hand, no nonlinear 
effect appears for a parallel beam. 

Of course, this nonlinearity is not a problem for the communication system: it can 
easily be seen that this vacuum effect is much smaller than that of the atmosphere. In fact, 
the non trivial question is the opposite: how can we detect such a nonlinearity? In order to 
have a significant effect, we need very high fields to compensate ^ = 6.7 x lO^^^m^/7. 
Which is the maximum electromagnetic energy density that can be achieved? 



ULTRAHIGH POWER LASER BEAMS 

Ultra intense photon sources are available thanks to the discovery of chirped pulse am- 
plification (CPA)LI6] in the late 80's and optical parametric chirped pulse amplification 
(OPCPA) ^ in the 90's. Very recently, the HERCULES Laser [0] reached a peak in- 



tensity / ~ 2 X lO^^Wcm^^, corresponding to a value i^p ~ 4 x 10^^^ for the product 
giving the importance of the nonlinear QED effects. In few years, the projected European 
Extreme Light Infrastructure (ELI) [8J will achieve / ~ lO^^Wcm^^, or i^p ~ 2 x 10^^. 
Finally, for a more distant future there is the hope to achieve / ~ 10^^ Wcm^^ [6], or 
~2x 10-^ 

These are the highest values for the intensities of the electromagnetic field that are 
available (HERCULES) or will be available in the future. These values are enormous. 
The power per pulse that has been reached at HERCULES is ^ 3 x lO^^W, which IS an 
order of magnitude larger than the total power used by Mankind at present! Of course, 
it is the ultrashort duration of the pulse (At ~ 3 x 10^ ^^5) that make it possible to reach 
such intensities keeping the energy per pulse as small as ~ 107. On the other hand, at 
ELI ^ the power per pulse will be lO^^VF, an order of magnitude larger then the total 
power that the earth receives from the sun! (In that case. At ~ 10^^^ — lO^^^i', so that 
the energy per pulse will still be reasonably small ~ 1 — 10-^7.) 



Of course, these facilities have very important technological and scientific applica- 
tions. First, it has been proved that they can be used to produce beams of high energy 
charged particles after bombarding a plasma [6]. This effect is highly nonlinear, since it 
can produce kinetic energies of the order of the MeVs or even the GeVs (and possibly 
even higher) using the ~ eV photons of the laser bumps This will provide cheap 
and reliable sources of radiation for cancer therapy and medical diagnostics, and may 
also help to ignite nuclear fusion. Moreover, it can be hoped that this surprising new 
form of accelerating particles can lead to a new era in particle physics. A second class of 
applications uses the small duration of the pulses, that will allow for an ultrashort time 
microscopy, which can be very important e.g. in nanotechnology. 

Finally, as we shall see below, these ultrahigh power lasers will provide a very 
promising tool for the exploration of the nonlinearity of the vacuum, and in particular 
for the search of 77 scattering. 



PREVIOUS PROPOSALS AND CURRENT LIMIT 

We have seen that the QED radiative corrections imply a nonlinear wave equations for 
the vector potential . On the other hand, the modified Maxwell's equations for the av- 



erage values of the electromagnetic quantum fields E and B were obtained long ago [lllll . 
and they were used to show that the QED vacuum exhibits the DC Kerr effect, implying 
birefringence fl^. More recently, several works have proposed different configurations 
that can actually be used for experimental tests of the nonlinear optical response of 
the vacuum, e.g. using harmonic generation in an inhomogeneous magnetic field [[l8|l . 
QED four- wave mixing[19], resonant interactions in microwave cavities [20], or laser- 
induced QED vacuum birefringence [21] which can be probed by x-ray pulses[22|]. 
among others [|23ll . All these proposals are very interesting and all of them deserve to 
be explored. They will provide different, independent tests for photon-photon scatter- 
ing. However, all these cases will need some technological advance. For instance, x-ray 
probing of QED birefringence requires an extra free electron laser, which is not yet 
available. Other techniques like four-wave mixing processes[19] require the crossing of 
at least three beams, with the corresponding alignment problems. As we shall see in the 
next section, our proposal is considerably simpler than the previous ones, and it only 
requires a ultrahigh energy pulse. 

Here, we would like to mention a very recent result obtained by the PVLAS collabora- 
tion [|24ll that has searched for an evidence of birefringence of the vacuum in a magnetic 
field background [25]. Their negative result has been used to set the current limit on 
the cross section for 77 scattering at optical wavelengths. Assuming the QED effective 
Lagrangian, but treating ^ as a free phenomenological parameter, their results can be 
translated in the limit i^exp < 3.1 x lO^^^m^/7. We see that this constraint is 4.6 x 10^ 
times higher than the QED value. Note that in Ref. [24] this result was also stated as a 
limit on the photon-photon cross section, Oyy < 4.6 x lO^^^m^, for A = 1064nm, and 
Oyy < 2.7 X IQ-^^m^, for X = 532nm. 



A NUMERICAL SOLUTION OF THE QED WAVE EQUATION 



In Ref. [14] we have obtained the first numerical solution of the full nonlinear wave 
equation Eq. (|5]) in the case of two counter-propagating plane waves that travel along the 
z-axis, for simplicity having the same phase at the space-time origin. The corresponding 
analytical solution of the linear wave equation (that can be obtained by setting i^/ = in 
Eq. ©) would be 

Mm{t,z) = Y[cos(fe-(DO+cos(fe + coO] (6) 
= £/cos((0?)cos(fe), 

where £/ is a constant amplitude, k = (0, 0, k) is the wave vector, and O) = ck is the 
angular frequency. It is easy to see that Eq. Q can also be considered as the analytical 
solution of the linear wave equation satisfying the boundary conditions 



A(?,0) = s^cos{cot), (7) 
<9,A(?,0) = 0, 

where for convenience we chose as the integration interval a 'small' cuboid of time- 
dimension 2n/(0 and space-dimension 2n/k. 




FIGURE 1. Numerical solution Anum/-i^ of Eqs. ^ and (jT) for = 0.00252, as a function of the 
adimensional time and space coordinates t = cot and C, =kz (from Ref. 13 )• 



In Ref. [|l4ll . we have found the solution of the nonlinear Eq. ([5]) that satisfies the same 
boundary conditions of Eq. (|7]), using the same small cuboid as the integration interval. 
This solution is shown in Figs. [T] and |2] for a choice of parameters such that (^£0^2/^(0^ = 
0.01. The corresponding time-averaged value p of the energy density turns out to be 

approximately constant along the z-evolution, giving E,p = -^t, j^^^^^pdt = 0.00252. 

As we have discussed above, this value of the product i^p is several order of magnitude 
larger than what can be achieved in the laboratory in the next decades. Of course, we 
will use realistic values of p when we will present our proposals of experiments in 



FIGURE 2. Zero-time plot of A num/jz/, as a function of the adimensional space coordinate ^ — kz (from 
Ref. Iil4i1). 



the last sections. For the moment, it is interesting to note that even such an enormous 
energy density is still small enough so that the effect of the nonlinear terms gives a small 
correction to the linear evolution in the short distance. In fact, we see from Fig. |3] that 
the relative difference between the linear and the nonlinear evolution, in our integration 
interval which is of the order of the wavelength, is of the order of few percent, i.e. of 
the same order than the adimensional parameter ^£o£/^(0^ = 0.01. This result is not 
surprising, and will provide a justification for the perturbatively-motivated variational 
approach that we will use in the next sections. 




FIGURE 3. Relative error (Anum ~^iin) of the linear approximation, as a function of the adimen- 
sional time and space coordinates T = cot and = kz (from Ref. (l^ ). 

However, from Fig. [3] itself, we can also appreciate that the difference between the 
linear and nonlinear behavior tends to increase along the z evolution, so that it can 
be expected that it will eventually become large after a distance much larger than the 
wavelength. In the next section, we will see an analytical argument that confirms this 
expectation. 

Finally, in Fig. |4] we compare the z-evolution of the solutions Anum and Aiin in a 
greater detail for values of z around the second zero of the solutions. We see that Anum 
anticipates Ajjn, and the corresponding phase shift can be evaluated numerically if we 
define an effective wave vector component by computing the value zo corresponding 
to Aiin(0,zo) = 0, and setting k^zo = 37r/2. The numerical determination of the zero 
gives ^0 = kzo = 4.68885, so that k^ = l.OOSOfc. We will provide a full explanation for 
this result in the following section. 
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FIGURE 4. Detail of the zero-time functions A-awal (upper curve) and Ann/^/ (lower curve), for 
values of the adimensional space coordinate ^ = kz close to the second zero (from Ref. |14]). 



Although it can be considered as an interesting achievement due to its simplicity and lack 
of previous approximations, in practice the numerical solution that has been discussed 
above can only be obtained in the special configuration of two counter-propagatin g, in - 
phase waves, and for short propagation (of the order of the wavelength). In Ref. PUD, 
we have found a variational approximation that can be used to provide an analytical 
solution which is valid even for a long evolution. In particular, we have considered the 
same two counterpropagating waves as in the previous section, that would be described 
by Eq. ^ in the linear case. Note that any of the two crossing waves, ^cos(fe— cot) 

and ^cos(fe-|- cot), when taken alone, would be a solution of both the linear and non- 
linear equations, provided that CO — ck. However, their superposition would only solve 
the linear equations of motion. Taking into account the smallness of the product ^p, and 
using the Variational Method, we have found the following analytical approximation 



where x = <^/£o-^^c^fc^/2 ~ 2^ipk is a parameter that describes the leading nonlinear 
effects. 

As a result, we find that after a distance Az the phase of the wave is shifted by a term 



Note also that the result of Eq. ([8]) can be stated equivalently by defining a wave vector 
as kz = k + X, that satisfies a modified dispersion relation, CO = c{k^ — x)- 

In QED {^i = for the same choice of parameters that was considered in the pre- 
vious section, — ^Cqs^^co^ /A = 0.0025, in agreement with the numerical simula- 
tion. In this case, we get x = 0.005^, which is two order of magnitude smaller than 
k. Our approximations can then be expected to be reasonably good even for the ex- 
tremely large of p that we have chosen here. Note also that this corresponds to a value 
kz = k + X = 1.005^, in agreement with the result that we obtained from the numerical 
simulation in the previous section. 

In Fig. [5] we compare the corresponding analytical solution Avar, given in Eq. ([8]), 
with the numerical solution Anum of the QED wave equations that we have found in 
the previous section. Comparing with Fig. [3l we see that the variational solution is an 



VARIATIONAL APPROXIMATION 



A = cos (cot) cos[{k + x)z], 



(8) 



= xAz ~ l^ipkAz 



(9) 



FIGURE 5. Relative error (Anum — ^var) /-^s^ of the variational approximation, as a function of the 
adimensional time and space coordinates T = cot and ^ = kz. 

order of magnitude closer to the numerical simulation than the linear evolution, Eq. 
that was obtained by completely neglecting the nonlinear terms. This is a significant 
improvement in such a short distance. However, according to the previous discussion, in 
the perturbative regime corresponding to small values of the product ^p, the variational 
solution is expected to be a good approximation even when a longer propagation distance 
is considered along the z-axis. 

On the other hand, the agreement of the variational solution with the numerical 
simulation can be used as an additional, a posteriori justification for our analytical 
approach. Therefore, in Ref. [14] we have applied the Variational Method to different 
configurations, that do not allow for a direct integration of the full Eq. The most 
promising one is the case of the scattering of a relatively low power wave with a 
counterpropagating high power wave, both traveling along the z axis. We have also 
allowed for an arbitrary initial phase difference between the two waves, assuming for 
simplicity that two waves have the same frequency. In this case, the main effect will be 
the modification of the low power beam due to the crossing with the high power one, 
which would be unaffected in a first approximation. Using our variational approach, we 
have found that the low power wave is modified, as compared to the linear behavior, as 
described by the following approximated analytical solution 

Ai{t,z) = aocos[{k + ri)z-cot], (10) 

where rj = 2^i£o£/^c^k^ ~ and p is the (average) energy density of the high 

power wave. Therefore, after a distance Az the low power beam accumulates a phase 
shift 

A4> = T]Az ~ 4(^;pA:Az. (11) 

We stress that this result only depends on the energy density of the high power wave, 
as far as it much larger than that of the low power wave. The initial phase difference is 
found to be irrelevant. 



PROPOSAL OF EXPERIMENTS AT HERCULES AND AT ELI 



In Ref. UM, we have proposed an experimental setup which can provide the simplest and 
most sensitive approach for the optical search of photon-photon scattering in vacuum. 
We have also showed that this porposal can be used at future exawatt facilities such as 
ELI ^ to detect the QED effect. Here, we will briefly review that idea, and show that it 
can also be used at present high power laser facilities, such as HERCULES jvl], to search 
for non-QED effects, or to improve the current limit for photon-photon scattering. 

The proposed experiment can be described as follows. A common laser pulse is 
divided in two beams, A and B, one of which (say A) crosses at a 180*^ a very high 
power beam. As a result, the central part of the distribution of the beam A has acquired 
a phase shift with respect to beam B, that has propagated freely. 

In an experiment corresponding to the parameters of the ELI project in its first step we 
have pulses of wavelength A = 800nm, intensity / = lO^^W cm^^ and duration A? = 10/^' 
which are focused in a spot of diameter d ~ 10/im. From Eq. ([TT]) . using for t,i the QED 
value of as given in Eq. dS]), this results in a phase shift A4> ^ 2 x lO^^raJ for beam A, 
which can be resolved comparing with the beam B which was not exposed to the effects 
of QED vacuum. Current techniques like spectrally resolved two-beam coupling, which 
can be applied for ultrashort pulses[26], can be used to this purpose, being sensible to 
phase shifts as small as A4> = 10^^ rad. It is interesting to note that the sensitivity of 
this method for the detection of photon-photon scattering may be enhanced by a suitable 
choice of the combination of the intensity /, the wavelength A and the time duration At 
that enter in Eq. (fTTI) . In fact, taking into account that Az — cAt, Eq. (fTTI) implies that the 
most sensitive experimental configuration will be that having the maximum value of the 
combination I At / A . 

Moreover, in an experiment corresponding to the parameters of the HERCULES laser 
[0], A = SlOnm, / = 2 X lO^^Wcm^^, T = 30 fs, d ^ 0.8/im, the same configuration 
as above would give A4> ~ 10^^ (from QED), which is possibly beyond the possibility 
of detection using current techniques. However, the experiment could find evidences 
of non-QED 77 scattering, or, if the result is negative, it can be used to substantially 
improve the current limit. Assuming as in Ref. ll24ll the same form of the effective 
Lagrangian than in QED, and taking A4> = lO^raJ as the sensitivity in the phase shift 
measurement, we see that HERCULES can prove values of the t,i parameter down to the 
value ~ 5 X 10^^^, which is almost two orders of magnitude lower than PVLAS current 
limit. Since the scattering cross section depends quadratically on t,, HERCULES can 
improve the current limit on by a factor (3 x 10^26/5 ^ 10-28)2 ^ 4 x 10^. This 
sensitivity may be increased by the upgrades of HERCULES that are programmed for 
the next months. Again, the most sensitive configuration will be that maximizing the 
product I At I A . 

Finally, it could be interesting to study the possibility of further enhancing the sensi- 
tivity of our proposed experiment by focalizing the high power pulse in the center of a 
cavity, keeping the intensity on the mirrors below the damage threshold. In this case, the 
high and low power pulse might be synchronized to cross each other several times. The 
number of times they meet would then multiply the value of the phase shift. 



CONCLUSIONS 



We have reviewed our recent proposal lfl3l [Till of a very simple experiment for the 
search of photon-photon scattering in vacuum. In this proposal, a low power laser pulse 
is made to cross with an ultrahigh power laser bump. 77 scattering will then be tested by 
measuring the phase shift of the low power laser beam, e.g. by comparing with a third 
low power laser beam. Even in the first step of ELI, we have found [ 1 4] that the resulting 
phase shift due to the QED-induced 77 scattering will be at least A^> 2 x lO^raJ, 
which can be measured with present technology. Moreover, in the present work we 
have also discussed the possibility of performing our proposed experiment at present 
ultrahigh power laser facilities, such as HERCULES. In particular, we have shown that 
HERCULES is already able either to detect 77 scattering, or to improve the current limit 
on the cross section Oyj at optical wavelengths by more than three orders of magnitude. 
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